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determined  by  neighboring  date.  Moreover,  the  effect  of  distant  data  on  algorithm  for  computing  plecevi.se  polynomial  re  present  all  ona  of  splines. 
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The  ter.  "local  dependence"  was  originally  applied  to  least-squares  In  ■"  of  thl»  «ork,  the  local  dependence  properties  of  spline 

spline  approximation  by  Powell  (P2)  In  a paper  bounding  the  effect  of  approximations  were  derived  from  the  corresponding  local  dependence 


lease-squares  splines  problems  1.  limited  storage. 


r 


This  dissertation  is  divided  into  six  chapters:  I)  Introduction 


pages  v and  vi 


Splines  and  B-Splines  a function  s(t)  in  the  collection  of  splines  S(k,tj  is  a linear 


to  uilty,  i.e.,  of  fc-l  order  B-spllnes,  where 


Increasing  k.  For  example,  with  infinitely  many  uniformly  spaced  knots 


the  discrete  Cram  matrix  la  singular,  and  the  discrete  THEOREM  4.2 


proof  la  somewhat  lengthy  [E5] 


■d  lately 


Proof ; Consider  the  continuous  real-valued  function 


COROLLARY  5.4 


otherwise  noted,  the  operations  counted  will  be  multiplications 


the  storage  counts 


TABLE  3«1  CUMMLNT  I Diagonal  elements  and  forward-solution 


» j:*mx(  i-pli)+pli-l J,  fl  ) UNTIL  i-1  DO  two-dioenaional  array  to  store  the  matrix,  because  many  FORTRAN  systems 

g[pf*)+i)  !-  g[pt*J+H  - g(pU1+J]  * g[J+p[tJ]  . , 

- require  an  integer  multiplication  for  each  two-dimensional  array 


ALCORITttt  4.3:  Local  Binary  Search  |B4J 


the  algorithm  reqvd  res  at  most  2N  log2(  K ) floating  point  comparisons. 


The  most  straightforward  evaluation  scheme  is  the  s u>  (II. 2. 14) 


spline  Is  a simple  matter  of  evaluating  the  k-*  order  spline  (5.3) 


S3UTT  dg-fl 


III. 6 The  Piecewise  Polynomial  Representation 


ALOORITW  6.1:  Coefficient  Piecewise  Polynomial  Evaluation  In  this  representation,  the  values  of  the  spline  and  its  derivatives  are 


very  efficient  incremental  evaluation  scheme  can  be  employed.  Suppose 


that  the  spline  Is  to  be  evaluated  at  the  points  «“1«  these  Iterates  are  generally  not  exact,  they  are  polynomials 

In  1,  l.e.,  (see  [Kl,»  l.  2. 6.  ( 9)  ,(40)1 ) 


evaluation  scheme  requires  k-I  additions  per  point.  for  al  1 but  certain  fortuitous  choices  of  initial  conditions 


Operation  counts  for  the  algorithms  described  in  this  section  are 


Algorithm  5.3.  The  second  algorithm  is  a table  look-up  algoritl 


[i  :«  i+1 ; 
u (nrt-1  J : » t (n+  1 } 


If  the  knots  are  uniformly  spaced 


ALGORITHM  7.3:  Conversion  to  Piecewise  Polynomial,  Translates  Operation  counts  lor  the  algorithms  described  In  this  section  are 


FIClffF.  B.2  F1CURE  8.3 

Operation  Counts  to  Evaluate  k ■ A B-splines  for  large  n an 
The  B-Spllnes  Evaluated  at  the  Knot  t^  with  (♦)  and  without  (■)  converting  to  piecewise  polynomials 


The  piecewise  polynomials  for  the  B-splines 
in  each  interval  through  the  PUT  function 


provided  chi*  series  converges.  However,  the  Neumann  series  (3.1)  Since  PA  P is  symmetric  and  positive-definite  (S«,p.  189, p.  30b] 
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COROLLARY  4.2  THEOREM 


Q 


COROLLARY  6.4  From  Theorem 


these  blocks  can  be  computed  using  a fixed  number  of  operations 


The  Cholesky  factor  L is  a lower  triangular , block  bl-dtagonal  and  from  (7.2)  we  can  obtain  an  expression  for  the  products  of  the 

matrix  with  the  same  partitioning  as  A.  Blocks  I.  t and  L.  ...  satisfy  diagonal  blocks  of  the  Cholesky  factor  L 


ll(Ali-l)'1)  a H »•  l±r+l 


problems  in  limited  storage.  Other  local  dependence  results,  though  not  [D3,D4]).  In  §3,  we  give  several  simple  numerical  examples  i'or 


Chapter  IV  to  develop  local  dependence  and  local  solution  results  for  which  the  Gram  matrix  has  an  arbitrarily  large  condition  number.  For 
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increases  with  |l“j|  and  sums  of  the  form  £ ®jl  <1°  not  converge. 


first  term  of  the  previous  result,  the  local  approx imat ion 
n be  bounded  using  Corollary  3.3.  If  r-i  is  sufficiently 


The  actual  damping  constant  2-^3  ■ . 26tt  is  considerably  smaller  than  our  Otherwise 


2 70.4  (.724) 


while  the  rth  row  of  L,  D,  and  £ is  being  computed,  the  arra>s 


results  are  stored  in  the  arrays  flk,k]  and  c[k).  For  example,  with 


estimates  for  £ can  be  obtained  from  (5.1)  using 


which  Is  equivalent  to  (for  e y > 0) 
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Nonzero  Entries  In  Rows  of  the  B-Spllne  Cram  MatrU  Th*  *2  Condltion  lju“ber  °l  the  B'*PliDe  trM‘  **lTlx 

for  Infinitely  Many  Uniformly  Spaced  Knots  for  Infinitely  Many  Uniformly  Spaced  Knots 


Alternative* — Par':  I.  International  Statistic*7  Review  42(1974), 
147-174. 


